Critical behavior at the interface between two systems belonging to different 

universality classes 



o 
o 

(N 
5h . 

< 
(N 



■4— > 

CO 

a 

C 

o 

> 

00 
(N 

o 

O 

c3 



13 

o 
o 



X 
S3 



Farkas A. Bagamery, 1,2 Loi'c Turban, 2 '* and Ferenc Igloi 3 ' 1 ^ 

1 Institute of Theoretical Physics, Szeged University, H-6720 Szeged, Hungary 
Laboratoire de Physique des Materiaux, Universite Henri Poincare (Nancy 1 ), BP 239, 
F-54506 Vandceuvre les Nancy Cedex, France 
3 Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O.Box 49, Hungary 

(Dated: February 6, 2008) 

We consider the critical behavior at an interface which separates two semi-infinite subsystems 
belonging to different universality classes, thus having different set of critical exponents, but having 
a common transition temperature. We solve this problem analytically in the frame of <j> k mean- 
field theory, which is then generalized using phenomenological scaling considerations. A large 
variety of interface critical behavior is obtained which is checked numerically on the example of 
two-dimensional g-state Potts models with 2 < q < 4. Weak interface couplings are generally ir- 
relevant, resulting in the same critical behavior at the interface as for a free surface. With strong 
interface couplings, the interface remains ordered at the bulk transition temperature. More inter- 
esting is the intermediate situation, the special interface transition, when the critical behavior at 
the interface involves new critical exponents, which however can be expressed in terms of the bulk 
and surface exponents of the two subsystems. We discuss also the smooth or discontinuous nature 
of the order parameter profile. 



I. INTRODUCTION 

Systems which undergo a second-order phase transi- 
tion display singularities in different physical observables 
which have been the subject of intensive research, both 
experimentally and theoretically. 1 At the critical tem- 
perature, T c , due to the existence of a diverging correla- 
tion length £ ~ \T — T c \~ u , microscopic inhomogeneities 
and single defects of finite size do not modify the critical 
singularities which are observed in the perfect systems. 2 
However, inhomogeneities of infinite extent, such as the 
surface of the sample, 3-5 internal defect planes, 6 etc., 
may modify the local critical properties near the inhomo- 
geneity, within a region with a characteristic size given 
by the correlation length. For example, the magnetiza- 
tion, m, which vanishes in the bulk as m ~ (T c — T)@ 
behaves as mi ~ (T c - Tf 1 at a free surface 3 " 5 and the 
two critical exponents, (3 and /3i, are generally different. 

Inhomogeneities having a more general form, such as 
localized 7 and extended defects, 8 corners, 9 wedges and 
edges, parabolic shapes, 10 etc., often have exotic local 
critical behavior; for a review, see Ref. 11. The local 
critical behavior can be nonuniversal, so that the local 
exponents vary continuously with some parameters, such 
as the opening angle of the corner, 9 the amplitude of a 
localized, 7 or extended defect. 8 The inhomogeneity can 
also reduce the local order to such an extent that the local 
magnetization vanishes with an essential singularity, as 
observed at the tip of a parabolic-shaped system. 10 On 
the contrary, for enhanced local couplings, a surface or 
an interface may remain ordered at or above the bulk 
critical temperature, 3-5 which in a two-dimensional (2D) 
system leads to a discontinuous local transition. 12 

In the problems we mentioned so far the inhomo- 
geneities are embedded into a critical system the bulk 
properties of which govern, among others, the divergence 



of the correlation length and the behavior of the order- 
parameter profile. There is, however, another class of 
problems, in which two (or more) systems meet at an in- 
terface, each having different type of bulk (and surface) 
critical properties. In this respect we can mention grain 
boundaries between two different materials or the inter- 
face between two immiscible liquids, etc. 

If the critical temperatures of the two subsystems are 
largely different, the nature of the transitions at the in- 
terface is expected to be the same as for a surface. 13 
At the lower critical temperature, due to the presence 
of the nearby ordered subsystem, the interface transi- 
tion has the same properties as the extraordinary sur- 
face transition. 3-5 At the upper critical temperature, the 
second subsystem being disordered, the interface transi- 
tion is actually an ordinary surface transition. 3-5 If the 
dimension of the system is larger than 2 and if the inter- 
face couplings are strong enough, one expects an interface 
transition in the presence of the two disordered subsys- 
tems whose properties should depend on the universality 
classes of these two subsystems. 

Even in 2D, the local critical behavior at the interface 
can be more complex if the critical temperatures of the 
subsystems are the same or if their difference is much 
smaller than the deviation from their mean value. In 
this case an interplay or competition between the two 
different bulk and surface critical behaviors can result in 
a completely new type of interface critical phenomena. 
In this paper we study this problem, assuming that the 
critical temperatures of the two subsystems are identical. 

The structure of the paper is the following. The mean- 
field solution of the problem including ip 3 , (p 4 and ip G 
theories and the interface between them is presented 
in Sec. II. The mean-field results are generalized in 
Sec. Ill using phenomenological scaling considerations. 
In Sec. IV these results are confronted with Monte Carlo 
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simulations in 2D for interfaces between subsystems be- 
longing to the universality classes of the Ising model, the 
three- and four-state Potts models, as well as the Baxter- 
Wu (BW) model. Our results are discussed in Sec. V and 
some details about the analytical mean-field calculations 
are given in the Appendixes A and B. 



II. MEAN-FIELD THEORY 
A. Properties of the ip k model 

1. Free energy 

We consider a system with volume V limited by a sur- 
face S in the Landau mean-field approximation. The 
total free energy is the sum of bulk and surface contribu- 
tions which are functional of the scalar order parameter 
(p(r) so that: 3-5 



(V) 



h M dV 



fs 



\dS. 



(1) 



(S) 



Near a second-order transition, the order parameter is 
small and the bulk free energy density fb[(fi] is written 
as an expansion in the order parameter and its gradient, 
limited to the following terms: 

f b [<p] = f b [0] + ^(V^) 2 + ^ 2 + h<p . (2) 

The second term, with C > 0, gives a positive contri- 
bution associated with the spatial variation of the order 
parameter. A = —at (a > 0, t — T c — T) is negative 
when T <T C and measures the deviation from the criti- 
cal point. The next term with B > ensures the stabil- 
ity of the system in the ordered phase. In the last term, 
h is the bulk external field. When k is odd, the order 
parameter is supposed to take only non-negative values; 
otherwise the system would be unstable. 

In the same way the surface free energy density is writ- 
ten as 



f s [if)] = f s [0] + 



9±^L 

2 A 



(3) 



where ip is the value of the order parameter on (S). The 
constant C' s is positive and A is a characteristic length 
related to the surface and bulk couplings of the corre- 
sponding microscopic Hamiltonian of the system. 3 



2. Ginzburg- Landau equation 

The mean-field equilibrium value of the order parame- 
ter, tp(r), minimizes the free energy in (1). It is obtained 
through a variational method by calculating (Ji* 1 [</?], the 
change of the free energy, which vanishes to first order 



in the deviation Sip(r) of the order parameter from its 
equilibrium value. Using Eqs. (l)-(3), one obtains 

6F[<p] = [ [CVif-VSip+ (Aip + Btp k - 1 -h)6tp]dV 

J(V) 

+ l s) ( Cs ^ S(PdS - W 

The first term in the volume integral may be rewritten 

as 

CVip ■ VSip = V • (CSpVp) - C<5^VV , (5) 

and the contribution to (4) of the first term on the right 
can be transformed into a surface integral through Gauss' 
theorem. Then 

5F[tp] = [ (-CV 2 <p + A(p + Btpt'- 1 -h)5<pdV 

J(V) 

+ [ (-Cn-V<p + C s ^)6cpdS, (6) 
J(S) v A/ 

where n is a unit vector normal to the surface and point- 
ing inside the system. 

At equilibrium, the first-order variation of the free 
energy vanishes. The volume integral leads to the 
Ginzburg-Landau equation 



-CV>(r) + Aip(r) + B^{r) = h(r 



(7) 



governing the equilibrium behavior of the order parame- 
ter in the volume of the system and the surface integral 
provides the boundary condition: 



Cn-V<p{r)\ {s) =C s ^ 



(8) 



3. Bulk critical behavior 

In the bulk, the first term in (7) vanishes. The zero- 
field magnetization vanishes when T > T c and is given 

by 



(9) 



in the ordered phase (T < T c , h = 0). 

The connected part of order-parameter two-point cor- 
relation function is given by 



(10) 



where ip(r) is the equilibrium order parameter, solution 
of the Ginzburg-Landau equation. Taking the functional 
derivative of Eq. (7), one obtains: 

-CV 2 r G(r,r') + [A+{k-i)B V k b - 2 ]G{ry) = k B TS(r-r') . 

(11) 
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This may be rewritten as 



ip A model, 1,3(0) < 1: 



(-Vl+C 2 )G(r, r i) 



k B T 
C 



S(r - r') . 



(12) 



where the expression of the bulk correlation length £ fol- 
lows from Eqs. (9) and (11) and reads 



i 

V= 2 



(13) 



in the ordered phase. 



4- Order parameter profiles 

We now assume that the surface of the system is lo- 
cated at z — so that ip = ip(z). Then, according to 
Eqs. (7), (9) and (13), the zero-field normalized order 
parameter profile <p> = ip/ifb is the solution of the follow- 
ing differential equation: 



^ A ,% J. B ,„fc-2,ftfc-l 
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k-l 



dz 2 



C 



(k-2)e 



(14) 



Multiplying by 2D(p/ dz and taking into account the bulk 
boundary condition, dip/dz — > when ip — > 1, a first 
integration leads to: 



dp 
dz 



20 k -kip 2 + k-2 
k{k-2)e 



(15) 



To go further we have to specify the value of k and to dis- 
tinguish between surfaces (or interfaces) which are more 
ordered (<p(0) > 1) or less ordered (ip(0) < 1) than the 
bulk. Below we list the solutions of Eq. (15) which will be 
needed in the sequel. We use the notation <p+(z) [<p_(z)] 
for a system located in the z > [z < 0] half-space. The 
values of the integration constants Z + and l_ are deter- 
mined by the boundary conditions at z = 0. 
<p 3 model, ip(0) > 1: 



<P±(z) 



d(p± 
dz 



3 cotli 
cosh 

2£ ± 



z±l ± 



sinh 



2U 



(16) 



p 3 model, ip(0) < 1: 



dz 



3 tanh 
±A sinh 



z±l ± 
2C ± 



- 1 



cosh 



1±. 

2£ ± 



(17) 



p A model, ip(0) > 1: 

d<p± 



z±, 



dz 



± coth 
T J-sinrT 



2£ ± 



2£± 



(18) 



ip±(z) — ±tanh 



z±L 



d(p± 
dz 



1 



±— cosh 



2£ ± 

2/ « 



2£, 



(19) 



p 6 model, ip(0) > 1: 



p ± {z) = V2 



3 tanh^ 



z±h 



1 



-1/2 



dz 



3/2 



<p 6 model, <£(0) < 1: 



(p±(z) — ±sinh 



dz 



i— COsh ( Jj- 



sinh' 
sinh 



2£ ± 

2 [ 

2£ ± 



(20) 



1-1/2 

-3/2 

•(21) 



B. Surface critical behavior 

In this section we briefly consider, for later use, the 
critical behavior at the surface when the bulk is in its 
ordered phase (T < T c ). We suppose that the system 
is located in the z > half-space. Since there is no 
ambiguity here we drop the index + so that, for example, 
p(z) stands for <p + (z). In this geometry, according to 
Eq. (8), the boundary condition reads 



<p(0) _ dip 

Cs HT- c Tz 



(22) 



Below we list the values obtained for the integration con- 
stant I and the surface order parameter ip(0). 



1. ip 3, surface 

a. A < 0. When A < 0, the surface remains ordered at 
the bulk critical point, which corresponds to the extraor- 
dinary surface transition. Thus, the order parameter pro- 
file is given by Eq. (16) and I is obtained by expanding 
both sides of Eq. (22) in powers of Z/£ <C 1. To leading 
order, one obtains 



and 



* = 2^|A| 



^° ) = 2{c\A\) ^ 



(23) 



(24) 



4 



b. A > 0. In this case the surface is less ordered than 
the bulk and we have an ordinary surface transition. The 
profile is given by Eq. (17) and the solution is obtained 
by assuming that the ratio Z/£ is a constant. Details of 
the calculation are given in Appendix A. The boundary 
condition in (22) is satisfied with 



/ = 2 tanh 



-1/2 



so that: 



¥>(0) 



I CA 

71 a& 



<p t 



3/2 



(25) 



(26) 



Thus the surface exponent, /3\ — 3/2, is larger than the 
bulk exponent, (3 = 1. 

c. A — > oo. The profile is given by either (16) or (17) 
with / — * oo. Then the order parameter is constant keeps 
its bulk value until the surface and <p(0) = <fb — Pot- 
We have a special surface transition, which corresponds 
to a multicritical point where the lines of ordinary and 
extraordinary transitions meet with the line of surface 
transition 3 in a (T, 1/A) diagram. 



3. ip 6 surface 

a. A < 0. Once more we have an extraordinary surface 
transition with a profile given by Eq. (20). As shown in 
Appendix A this is another instance where the boundary 
condition in (22) is satisfied with a constant value of the 
ratio Thus, as in (25), we have 



I = 2 tanh" 1 ^-^=) £o^ 1/2 



(31) 



The leading contribution to the surface order parameter 
is then 



(32) 



b. A > 0. The profile at the ordinary surface transition 
is given by Eq. (21). Here we have the standard behavior, 
l/£ < 1, with 



I 



C_ 



A. 



(33) 



2. (p 4 surface 



The surface order parameter displays the following be- 
havior: 



a. A < 0. This corresponds as above to the extraordi- 
nary transition where the surface remains ordered at the 
bulk critical point. The profile is given by Eq. (18) and 
the boundary condition in (22) requires i/( < 1 so that 
one obtains 



'=!' A >' 



(27) 



The leading contribution to the surface order parameter 
is given by 



(28) 



b. A > 0. At the ordinary surface transition, the 
profile is given by Eq. (19). Here the boundary condition 
is satisfied with l/£ <C 1, which gives 



The surface order parameter vanishes as 

, . 1 CA 
2 C s c,o 



(29) 



(30) 



Thus, the surface exponent is f3\ = 1 to be compared to 
the bulk exponent, [3 = 1/2. 

c. A — > oo. Here too, the boundary condition in (22) 
leads to I — > oo and the surface order parameter keeps 
the bulk value, ip(0) = (prf 1 / 2 , at the special surface tran- 
sition. 



1 CA 

V6 Cs^o 



9o 



t 3 / 4 . 



(34) 



Thus the surface exponent is (3\ = 3/4 whereas [3=1/4 
in the bulk, 
c. A — > oo. 

As for the other models, the length I is infinite at the 
special transition and the surface order parameter has 
the bulk value, tp(0) = (pot 1 / 4 . 

Although the characteristic length / sometimes remains 
finite and sometimes diverges at the critical point, the 
exponent (3\ at the ordinary surface transition always 
satisfies the scaling relation: 

(3\= (3 + v (ordinary transition) . (35) 

In the same way, at the special transition, we have: 



(3\= (3 (special transition) . 



(36) 



One should notice that these scaling relations arc only 
valid in mean-field theory. 3 

At the extraordinary transition the singular term, gov- 
erning the approach to the constant value at T c of the 
surface order parameter, appears at the next order in 
the expansion. It vanishes linearly in t for the if 3 and 
tp 4 models and as t 1 / 2 for the <p e model. We do not give 
further details since we shall not need it in the following. 
For the same reason, we did not examine the properties of 
the surface transition which occurs in the surface region, 
above the bulk critical temperature, when A < 0. 
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Interface critical behavior 



have 



In this section, we consider the critical behavior at 
the interface between two systems, belonging to differ- 
ent universality classes, in their ordered phase (T < T c ). 
Thus the free energy densities of the two subsystems are 
given by (2) with different values of k. They are coupled 
through an interface at z — with free energy density 



(37) 



where 



I- = /§-|A|, 



(39) 



A* 

/=1+ \ 1+3 ^' 



A* = C + Clip o fe>) 2 (aq) 



We assume that the positive half-space corresponds to 
the system which is the more ordered in the bulk when 
T c is approached from below so that /3 + < /3~ . The order 
parameter profiles, <p+(z) for z > and <f-(z) for z < 0, 
have now to satisfy 



The leading contribution to the order parameter at the 
interface, 



(41) 



V(0) 

y(Q) 

A 



<p + (0) = <f-(0).. 



dz 



C_ 



dz 



(38) 



is also independent of t; i.e., the interface remains ordered 
at the bulk critical point. 

According to (40) and (41), the asymptotic dependence 
on |A| is the following: 



These boundary conditions generalize Eq. (22) for the 
interface geometry where each subsystem contributes a 
normal derivative to the surface integral in Eq. (6). 

When two subsystems are coupled, the boundary con- 
ditions in Eq. (38) determine the integration constants 
l± — i.e., the complete order parameter profile. In the 
following, we give these integration constants as well as 
f(0), the value of order parameter at the interface, for the 
different types of interface considered. Technical details 
about the calculations can be found in Appendix B. 

As in the surface case, depending on the value of A, 
different types of interface critical behaviors are obtained 
(see Figs. 1-3). When A < 0, the interface remains or- 
dered at the bulk critical point and we have an extraor- 
dinary interface transition. When d > 2, the local or- 
der persists above the bulk critical temperature until a 
A-dependant interface transition temperature is reached. 
This transition, which always occurs in mean-field theory, 
will not be discussed further here. When A > the inter- 
face order parameter vanishes at the bulk T c as a power of 
t. This corresponds to the interface ordinary transition. 
When parameterized by 1/A, these two transition lines 
meet, together with the interface transition line when it 
exists, at a multicritical point corresponding to the spe- 
cial interface transition located at 1/A = 0, T = T c . 



tp(0) oc 



|A|- 2 |A|»A* 
jAj- 1 |A| < A* 



(42) 



b. A > 0. This corresponds to weak couplings between 
the two subsystems. When 



0<A<A C = 2^%+, 



(43) 



the order parameter decreases from both sides towards 
the interface. Then <f-(z) is given by (17) and <p+(z) 
by (19) with 



(44) 



l_ diverges when A = A c . Then f-(z) is a constant 
and keeps its bulk value for any z < 0. 

When A c < A < 00, f + (z) is still given by (19) and l + 
keeps the value given in Eq. (44) but now (fl-{z), which 
increases, is given by Eq. (16) with 



L=2coth- ,i±^ Wr" 



(45) 



ip 4 interface 



a. A < 0. This corresponds to strong couplings at 
the interface. The order parameter increases when the 
interface is approached so that <f-(z) and <p+(z) are given 
by (16) and (18), respectively. To leading order in t, we 



For < A < 00, the order parameter at the interface is 
always given by 



1 C+A , 



(46) 



6 



0.05 




0.10 




FIG. 1: Reduced order parameter profile (j) = f{ z )/Vb at 
the <p 3 — (p 4 interface as a function of £ = z/£+ for different 
values of A = A/£+ with 1/A = -0.2, -0.1, 0, 0.1 and 0.5, 
from top to bottom. The behavior for A > and £ < is 
enlarged in the inset where 1/A = 0, 0.1, 0.5, 1 and 10, from 
top to bottom. One may notice the change of behavior at 
1/A C = £d/A c = 1/2. All other parameters have the same 



values in the two subsystems: C + 
$ = £o and t = 10~ 4 . 



<Pt) : 



FIG. 2: Reduced order parameter profile <j> — ip(z)/ip£ at the 
ip 4 — ip e> interface as a function of £ = z/£ + for different values 
of A = A/£+ with 1/A = -0.1, -0.05, -0.03, and 0.1, from 
top to bottom. The behavior for A > and £ < is enlarged 
in the inset where 1/A = 0, 0.01, 0.02 and 0.1, from top to 
bottom. All other parameters have the same values in the 
two subsystems and t = 10~ 6 . 



where 



c. A — > oo. The profile remains monotonously increas- 
ing and keeps the same functional form as for A > A c 
although l_ and l + are now given by 



'- " 2 ( 3 tf«f 3(£ ° )2/3< 



so that 



2/3^-1/3 
2/3 



-1/3 



(47) 



2/3 



2. ip 4 — ip 6 interface 

a. A < 0. The interface is more ordered than the bulk. 
Thus <P-(z) is given by (18) and <p+(z) by (20) with 



I. 



l + = 2 



tanh 



V3 



+ h(t) 



+,-1/2 



e t 



32 V Ci^o Co / 



(49) 




A* 



(50) 



Here and below in Eq. (57) we keep the next-to-leading 
term h(t) in l + . This correction is actually needed to 
obtain the correct form of the profile in the vicinity of 
z = 0. The leading contribution to the order parameter 
at the interface is constant: 



fC_\A\ 



t y\ (48) 

Its asymptotic behavior 



(p(0) <x 



|A|" 
|A|" 



1 |A| » A* 
V 2 |A| < A* 



(51) 



(52) 



follows from Eqs. (50) and (51) 

b. A > 0. The profile is always decreasing when the 
interface is approached. Thus <P-{z) is given by (19) 
and <p+(z) by (21) with the following expressions for the 
integration constants: 



I- = \ t 



2 C + y^o 

3 c^ e + 



At 



-1/4 



Ci 



A. 



(53) 
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The interface order parameter behaves as 

1 C + A 



0.05 



<p(0) 



(54) 



c. A — > oo. 

Then the profile increases monotonously with z. ip_ (z) 
is given by (18) and <p+(z) by (21) with 



l_ = \2V6 



CM 



+ €0 £d 



1/2 



1/2 



t 



-3/8 



-3/8 



The interface order parameter vanishes as 

\ 1/2 



3 c_a 



(55) 



(56) 




5. <^ 3 — ip 6 interface 

a. A < 0. As usual in this case, the interface is more 
ordered than the bulk. The profiles, <fi-(z) and <p+(z) 
are given by Eqs. (16) and (20). The calculation of l_ 
involves the solution of an equation of the fourth degree 
(see Appendix B). Here we only report the limiting be- 
havior for large and small values of |A|: 



2 i~ |A| ' 



|A| > A*, 



L - (2V27 



;C+|A| \ ip Q ^ 



^ Co, |A|«A% 



h(t) 
h(t) 



= 2 



tanh 1 



V3 



h(t) 



-1/2 



2V3 /C-|A|V 



+ \ 2 
|A| « A* . 



|A| > A* 



4 aa 

The crossover is taking place around 



A* = C, 



2/3 x 1/3 

fo_ \ ( c +\ (to 



The interface order parameter reads 

2 



4/3 



(57) 



(58) 



<P{0) 
<p(0) 



3 / Cigp 
2 I C*_A 



2V3 



C + |A| 



|A| » A* , 
|A| « A* 



(59) 



b. A > 0. The profile is monotonously increas- 
ing. ifi-(z) and <p+(z) have the form given in Eqs. (16) 



FIG. 3: Reduced order parameter profile <j> — <p(z)/<p~£ at the 
ip A — ip e> interface as a function of £ = z/£+ for different values 
of A = A/£+ with 1/A = -0.2, -0.05, -0.02, and 1, from 
top to bottom. The behavior for A > and £ < is enlarged 
in the inset where 1/A = 0, 0.2, 0.5, 1 and 10, from top to 
bottom. All other parameters have the same values in the 
two subsystems and t = 10~ 3 . 



and (21) with the following values of the constants: 



L 



6V6. 



Co"*" 378 . 



The interface order parameter vanishes as 

1 C+A 



V(0) 



-^ 3/4 - 



(60) 



(61) 



c. A — > 00. The profile is still given by Eqs. (16) 
and (21) with the following values of the constants: 

1/3 



i2V6&4(c -) 2 a 



CM 



^oKo 



1/3 



-1/4 



At the interface we obtain 
tp(0) 



2 cm 



2/3 



Mo) 1/3 t 1/2 - 



(62) 



(63) 



III. SCALING CONSIDERATIONS 

Here we generalize the mean-field results obtained 
in the previous section. First, we consider the order- 
parameter profiles in semi-infinite systems with free and 
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fixed boundary conditions. These results are used after- 
wards to study the scaling behavior at an interface, which 
separates two different semi-infinite systems. 

A. Order-parameter profiles in semi-infinite 

systems 

We consider a semi-infinite system, which is located in 
the half-space z > 0, and which is in its bulk-ordered 
phase (T < T c ); see in Sec. II B. As in mean-field theory, 
the order-parameter ip(z) depends on the distance from 
the surface, z, and approaches its bulk value, (fib ~ 
for z/% ~> 1. The bulk correlation length asymptotically 
behaves as £ ~ \t\~ v . These expressions generalize the 
mean-field results in Eqs. (9) and (13). 

1. Free boundary conditions 

At a free surface, due to the missing bonds the local 
order is weaker than in the bulk. The surface order pa- 
rameter displays the so-called ordinary transition with 
the temperature dependence <p(0) ~ t" 1 , where generally 
[3i > (3. The profile, <p(z), which interpolates between 
the surface and the bulk value has the scaling form 3-5 

<p(z) = (fibford ' ( 64 ) 

and the scaling function, / or d(j/)> behaves as y^~^l v , 
for y <C 1. 



2. Fixed boundary conditions 

For fixed boundary conditions, the system displays the 
extraordinary surface transition and stays ordered in the 
surface region at the bulk critical temperature, so that 
ip(z) = 0(1) as t — > + and z<(. This behavior is for- 
mally equivalent to having a surface exponent, /3i = 0. 
The magnetization profile can be written into an analo- 
gous form 3-5 as in Eq. (64): 

<p(z) = fbfcxt ' ( 65 ) 

however, now the scaling function, f cy ±(y) has the asymp- 
totic behavior, 14 f e xt(y) ~ y~^^ v , for y <C 1. 

B. Interface critical behavior 

Now we join the two semi-infinite systems and study 
the behavior of the order-parameter in the vicinity of 
the interface. In general we expect that, depending on 
the strength of the interface coupling, at the bulk critical 
temperature the interface (i) can stay disordered for weak 



couplings, which corresponds to the A > case in mean- 
field theory or (ii) can stay ordered for stronger couplings, 
which is the case for A < in mean-field theory. These 
two regimes of interface criticality are expected to be sep- 
arated by a special transition point, which corresponds 
to A — > oo in mean-field theory 

To construct the order-parameter profile we start with 
the profiles in the semi-infinite systems and join them. 
First we require continuous behavior of the profile at 
z = 0, like in mean-field theory. The second condition in 
mean-field theory in Eq. (38) cannot be directly trans- 
lated; here we just use its consequencies for the extrapo- 
lation lengths. In the weak- and strong-coupling regimes 
in mean-field theory the left-hand side of Eq. (38) is fi- 
nite so that the derivative of the profile is discontinuous 
at z = and at least one of the extrapolation lengths 
l± is 0(1). The same behavior of l ± is expected to hold 
in scaling theory, too. On the other hand, at the spe- 
cial transition point in mean-field theory the left-hand 
side of Eq. (38) is zero and the extrapolation lengths are 
divergent. In scaling theory the asymptotic form of the 
extrapolation lengths is expected to be deduced from the 
same condition — - i.e. from the equality of the derivatives 
of the profiles. This leads to a relation l± <~ t~ v±UJ± , in 
which < uj± < 1 is defined later. 

If the subsystem — say at z > — has an ordinary tran- 
sition the interface magnetization follows from Eq. (64) 
as 

^(0)~i ft , i = (l-w + )ft+w + l3 + , (66) 

and (3 + < Pi < j3^. On the other hand, if the 
subsystem — say at z < — has an extraordinary tran- 
sition the interface magnetization exponent follows from 
Eq. (65) as 

Pi =lj_(3_. (67) 

Evidently, [3i calculated form the two joined subsystems 
should have the same value. This type of construction of 
the order-parameter profiles will lead to a smooth profile 
at the interface provided the extrapolation lengths are 
smaller or, at most, of the same order than the correlation 
lengths, max(Z + ,Z_) < min(£; + ,£_), which holds provided 

max(w + f + ,w_i/_) < min(v +1 v_) . (68) 

Otherwise, the profile measured in a length scale, 
min(£; + ,£_), has a sharp variation at the interface and 
as the critical temperature is approached the profile be- 
comes discontinuous. Note that in mean-field theory, 
with v + — v_ and uj± < 1, the profile is always smooth. 

1. Relevance-irrelevance criterion 

Here we generalize the relevance-irrelevance criterion 
known to hold at an internal defect plane with weak 
defect couplings. 6 If two different critical systems are 
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weakly coupled, the operator corresponding to the junc- 
tion is the product of the two surface magnetization oper- 
ators. Consequently, its anomalous dimension Xi is given 
by the sum of the dimensions of the two surface opera- 
whcrc xf = pf /v±. Then the scaling 
j/i, in a (i-dimensional system is 



tors, Xi — x\ + x{ 
exponent of the defect, 
given by 



Vi = d t 



d-l 



Pt 0i 



(69) 



where di = d — 1 is the dimension of the interface. 

For yi < the weak interface coupling is irrelevant 
so that the defect coupling rcnormalizes to zero and the 
defect acts as a cut in the system. Consequently the in- 
terface critical behavior is the same as in the uncoupled 
semi-infinite systems and the interface magnetization ex- 
ponent is Pi — mm(0f) since the stronger local order 
manifests itself at the interface. In the other case, yi > 0, 
the coupling at the interface is relevant and the interface 
critical behavior is expected to be controlled by a new 
fixed point. 

For the 2D g-state Potts model with 2 < q < 4 we 
have 15 xi > 1/2; thus weak interface coupling is expected 
to be irrelevant according to Eq. (69). 

In mean-field theory, when d appears in a scaling rela- 
tion, it has to be replaced by the upper critical dimension 
d c for which hyperscaling is verified. However we have 
here different values of d c for the two subsystems so that 
there is some ambiguity for the value of d in Eq. (69). 
The analytical results of Sec. II C show that a weak in- 
terface coupling is also irrelevant in all the cases studied 
in mean-field theory. 



2. Weakly coupled systems 

For weak interface coupling the order parameter pro- 
file is not expected to display a maximum at the inter- 
face. Depending on the relative values of the critical 
exponents, P + , Pt, Pi it can be either a minimum, 
or an intermediate point of a monotonously increasing 
profile. We use the same convention as in Sec. II C, that 
P + < p~ and treat separately the different cases. 

a. [3 + < (3_ < Pt < Pi . The order-parameter profile 
is obtained by joining two ordinary surface profiles in 
Eq. (64) both for z < and for z > 0. In this case the 
weak coupling does not modify the asymptotic behavior 
of the more ordered, z > subsystem. Consequently we 
have ip(0) ~ (pt, l + = O(l) and u + = 0; thus Pi = Pt- 
From Eq. (66) we obtain 



U_ = 



Pi ~ Pt 
Pi - P_ 



(70) 



Note that the above reasoning leads to a smooth order- 
parameter profile at the interface, if according to Eq. (68) 
we have w_ v_ < i> + . This type of behavior is realized in 



mean-field theory for the ip 4 — ip 6 interface for A > 0, see 
in Sec. IIC 2. 

b. P + < P_ < Pi < Pi . In this case the profile is ob- 
tained from two ordinary subprofiles. The order parame- 
ter is still minimum at the interface, but it is determined 
by the z < subsystem, which has the larger surface or- 
der parameter. Consequently, tp(0) ~ <pl, l_ — 0(1) and 
cj_ = 0; thus Pi — Pl . From Eq. (66) we obtain 



Pt ~ Pi 
Pt-P + ' 



(71) 



and the order-parameter profile is smooth if ulj + v + < v_. 
This type of behavior is never realized in mean-field the- 
ory; see the exponent relation in Eq. (35). 

c. p + < Pi < P- < Pi- In this case the order- 
parameter profile is monotonously increasing and ob- 
tained by joining an extraordinary profile in Eq. (65) for 
z < with an ordinary profile in Eq. (64) for z > 0. The 
order parameter at the interface is determined by the 
surface order parameter of the z > subsystem. Then 
we have <p(0) ~ tpf, l + — 0(1) and u + = 0; thus Pi — p^. 
From Eq. (66) we obtain uo_ = P\jP- and the interface is 
smooth, provided < v + . In mean-field theory this 

type of behavior is realized for the ip 3 — ip 6 interface for 
A > 0; see in Sec. IIC 3. 



3. Special transition point 

In this case the profile is monotonously increasing and 
it is constructed by joining an extraordinary subprofile in 
Eq. (65) for z < with an ordinary subprofile in Eq. (64) 
for z > 0. As we argued before, the extrapolation lengths 
and the corresponding exponents are obtained (i) from 
the continuity of the profile at z = 0, 

ft = P-U- = (1 - u + )Pt + P + u + , (72) 

and (ii) from the continuity of the derivative at z = 0, 
which leads to the condition l + ~ l_ — consequently, 

Vi = = u)_v_ . (73) 

The solution of Eqs. (72) and (73) is given by 



Pr 



0- 



(74) 



Let us now analyze the condition for the smooth or dis- 
continuous nature of the interface given in Eq. (68). 

a. v_ > v + . In this case the condition is equivalent to 
P-/v- > P+/v + . As we will discuss in Sec. IV this con- 
dition is satisfied in 2D for the three- and the four-state 
Potts (or BW) models so that the profile is predicted to 
be smooth. On the contrary for the Ising and the three- 
state Potts models this condition does not hold; thus, the 
profile is probably sharp and becomes discontinuous at 
the critical temperature. Finally, for the Ising and the 
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BW models the relation in Eq. (68) is just an equality, 
so that we are in a marginal situation. 

b. v_ < v + . In this case the profile is smooth, provided 



1.00 



ft+<^T 



1 



(75) 



This type of situation seems to be less common in real 
systems. 



4- Strongly coupled systems 

In this case the interface stays ordered at the bulk crit- 
ical temperature, so that the profile is expected to be 
composed from two extraordinary subprofilcs. As a con- 
sequence the interface critical behavior is the same as in 
two independent semi-infinite systems, both having an 
extraordinary surface transition. 



IV. NUMERICAL INVESTIGATIONS 

We have studied numerically the critical behavior at 
the interface between two g-state Potts models on the 
square lattice, with different values of q for the two sub- 
systems. For a review of the Potts model, see Rcf. 16. 
In particular we considered the value q — 2, which corre- 
sponds the Ising model, as well as q — 3 and q = 4. All 
these systems display a second-order phase transition for 
a value of the coupling given by e qKc = 1 + ^/q, which 
follows from self-duality. The associated critical expo- 
nents are exactly known 18 for q = 2 (/3 = 1/8, v = 1 and 
ft = 1/2) and has been conjectured for q = 3 (/? = 1/9, 
v = 5/6 and ft = 5/9) and q = 4 ((3 = 1/12, u = 2/3 and 
ft = 2/3), where they follow from conformal invariance 20 
and the Coulomb-gas mapping. 21 We have also consid- 
ered the BW model, 17 which is a triangular lattice Ising 
model with three-spin interactions on all the triangular 
faces. This model is also self-dual and has the same crit- 
ical coupling as the Ising model. This model is exactly 
solved, 18 and it belongs to the universality class of the 
q = 4 state Potts model, but without logarithmic cor- 
rections to scaling, which facilitates the analysis of the 
numerical data. 

We have performed Monte Carlo simulations on 2D 
systems consisting of two LxL subsystems which interact 
directly through interface couplings Ki (between adja- 
cent spins of the two subsystems) such that Ki/K + = A. 
Here K + is the coupling in the half-space z > 0, which 
corresponds to the subsystem having the larger value 
of q, thus the larger magnetization. Periodic bound- 
ary conditions are applied in both directions. Using the 
Swendsen-Wang cluster-flip algorithm 19 we have calcu- 
lated the magnetization profile in systems with size up 
to L = 300 for different values of the reduced tempera- 
ture t = (T c — T)/T c and coupling ratio A. Depending 
on the size of the system and the temperature we have 
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FIG. 4: Order parameter profiles of the q = 3 — 4 interface at 
the reduced temperature t = 0.06 and for different values of 
the interface coupling ratio, A = 2.0,l.l,1.0, 0.9 and 0.5 from 
top to bottom. Inset: profiles around the special transition 
point, A = 1, for different temperatures. The results indicate 
a smooth profile at the transition point. 



skipped the first 5 — 20 x 10 4 thermalization steps and the 
thermal averages were taken over 6 — 20 x 10 6 MC steps. 
We have checked that the magnetization profiles, at the 
reduced temperatures we used, does not show any notice- 
able finite-size effects. From the magnetization data at 
the interface we have calculated effective, temperature- 
dependent interface exponents given by 



ft(i) 



log[m(t + S)/m(t - 5)} 
log[(t + S)/(t-S)} 



(76) 



which approach the true exponents as 5 —* and t — > 0. 
For the Ising-BW interface, we have also made calcula- 
tions at the critical temperature in order to check the 
finite-size scaling properties of the profiles. In the follow- 
ing we present the numerical results for the q = 2 — 3, 
q = 3 — 4 and Ising-BW interfaces. In each case we have 
a different type of special transition, separating the ordi- 
nary and the extraordinary transition regimes. 



A. q = 3 — 4 interface 

We start in Fig. 4 with a presentation of the order- 
parameter profiles, in the vicinity of the critical tem- 
perature, for different values of the interface coupling. 
Here one can differentiate between the ordinary transi- 
tion regime for small A, in which the magnetization at 
the interface vanishes faster than in the bulk of the two 
subsystems, and the extraordinary transition regime for 
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FIG. 5: Effective magnetization exponents measured on the 
two sides of the q = 3 — 4 interface for three different values of 
the coupling ratio, A. At the special transition point, A ~ 1, 
the theoretical prediction from Eq. (74) is indicated by a bar. 



large A, where the interface magnetization keeps a fi- 
nite value. The special transition separating these two 
regimes is located at A rj 1 The inset of Fig. 4 shows 
the evolution of the interface at the special transition 
point as the bulk transition point is approached. Here 
the criterion in Eq. (68) is satisfied, since, as discussed 
below Eq. (74), (3_/v_ = 2/15 > (3+/v + = 1/8. Thus the 
profile is predicted to be smooth, which is in accordance 
with the numerical results. 

The values of the effective, temperature-dependent ex- 
ponents, as defined in Eq. (76), are presented in Fig. 5 
for three values of A, corresponding to the different tran- 
sition regimes. Clearly the values of the effective expo- 
nents are affected by strong crossover effects for small t, 
since the limiting values are = (3\{q = 3) = 5/9 for 
A < 1 and = for A > 1, according to scaling the- 
ory. Unfortunately, due to finite-size effects we could not 
go closer to the critical point. At the special transition 
point, however, the crossover effects are weaker and the 
effective exponents are close to the theoretical prediction 
in Eq. (74), & = 32/363 = 0.088. 



B. q = 1 — 3 interface 

We have performed a similar investigation for the in- 
terface critical behavior of the q = 2 — 3 system and the 
results are summarized in Figs. 6 and 7. Here one can also 
identify the ordinary and the extraordinary transition 
regimes (see Fig. 6), which are separated by the special 
transition around A w 0.85. However, as can be seen in 
the inset of Fig. 6, the behavior around the special tran- 



FIG. 6: The same as in Fig. 4 for the q — 2 — 3 interface with 
interface coupling ratio A = 2.0, 1.0, 0.9, 0.85, 0.8 and 0.5 from 
top to bottom and t = 0.06. Inset: profiles around the special 
transition point, A c = 0.85, for different temperatures. The 
results indicate that the profile becomes discontinuous at the 
critical temperature. 
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FIG. 7: The same as in Fig. 5 for the q = 2 — 3 interface. 



sition point is more complex than for the q = 2 — 3 inter- 
face. The evolution of the profile suggests the existence of 
a discontinuity at the transition temperature. This is in 
accordance with the scaling criterion in Eq. (68) since, as 
discussed below Eq. (74), f3_/v_ = 1/8 < = 2/15 

leads to a discontinuous profile. Due to this disconti- 
nuity, it is more difficult to locate precisely the special 
transition point and to determine the associated interface 
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FIG. 8: Critical magnetization profiles in the Ising-BW sys- 
tem with two symmetrically placed interfaces for different val- 
ues of the coupling ratio A. 

exponent The measured effective interface exponents 
are shown in Fig. 7 for three values of A corresponding 
to the different interface fixed points. The crossover ef- 
fects are strong but, at the special transition point, our 
estimates are compatible with the scaling prediction in 
Eq. (74), fa = 25/237 = 0.105. 

C. Ising-BW interface 

The interface between the Ising model and the BW 
model (or the four-state Potts model) has some special 
features. These are mainly due to the fact that the 
anomalous dimension of the bulk magnetization in the 
two systems has the same value, f3-jv_ — (3 + /v + = 1/8. 
Consequently one can define and numerically study the 
finite-size scaling properties of the magnetization profile 
at the phase-transition point, since it is expected to scale 
as m{z,L) — L _1 / 8 /(^/L). The scaling function f(y) is 
expected to depend on the value of the interface coupling 
ratio A and we have studied this quantity numerically. 

The magnetization profiles at the critical temperature 
for different values of the interface coupling ratio A are 
given in Fig. 8. It is interesting to notice that the shape 
of the curves as well as the relative heights of the profiles 
in the two subsystems vary with the interface coupling. 
For A < A c w 1, the interface stays disordered and the 
interface critical behavior is governed by the surface ex- 
ponent of the Ising model. The larger bulk value on the 
Ising side is understandable since the profile on the right 
side is more singular, /3{" jv_ > (3^/u + ; see below Eq.(64). 
On the contrary, for A > A c the interface is ordered at 
the bulk critical temperature and the profiles decay to- 



FIG. 9: Scaled magnetization profiles for the coupled Ising- 
BW systems at the common critical temperature for L = 90 
- 300. The interface coupling ratio is fixed at the critical 
value, A C (L), for which the two maxima of the curves are 
identical. The inset gives the effective critical interface cou- 
pling as a function of the inverse size. 

wards the bulk values. 

At the special transition point A = A c , the profile 
has a universal form in terms of L 1 / 8 m(z/L). This is 
illustrated in Fig. 9, in which for each finite system a 
critical value A C (L) is calculated from the condition that 
the two maxima of the curves have identical values and 
the profile is measured at that interface coupling. The 
size-dependent effective interface coupling ratio, A C (L), 
shown in the inset of Fig. 9, seems to tend to a limiting 
value, A c w 1. The scaled magnetization profiles have 
different characteristics in the two subsystems. In the 
BW model, having the smaller correlation length, the 
profile has a smooth variation. On the contrary, on the 
Ising side, the profile has a quasidiscontinuous nature at 
the interface, which is probably related to the fact that 
in the criterion of Eq. (68) the equality holds. 

V. DISCUSSION 

In this paper we have studied the critical behavior 
at the interface between two subsystems displaying a 
second-order phase transition. We assumed that the crit- 
ical temperatures are identical but the sets of critical ex- 
ponents (i.e., the universality classes of the transitions) 
arc different for the two subsystems. By varying the 
interface couplings, we monitored the order at the in- 
terface and studied the behavior of the order-parameter 
profile as the critical temperature is approached. We 
provided a detailed analytical solution of the problem 
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in the framework of mean-field theory, which leads to a 
physical picture which is useful for the study of realistic 
systems. Solutions of the mean-field equations are ob- 
tained by adjusting the order-parameter profiles of the 
two semi-infinite subsystems through the introduction of 
appropriate extrapolation lengths on the two sides. The 
same strategy has been applied in the frame of a phe- 
nomenological scaling approach. As a result, basically 
three types of interface critical behavior are observed. 
For weak interface couplings the interface renormalizcs 
to an effective cut and we are left with the surface criti- 
cal behavior of the subsystems. In the limit of strong in- 
terface couplings, the renormalization leads to infinitely 
strong local couplings and thus interface order at the bulk 
critical point. Finally, for some intermediate value of the 
interface couplings, the interface displays a special tran- 
sition, which is characterized by a new critical exponent 
for the order parameter in the interface region. In the 
scaling theory this exponent can be expressed in terms 
of the bulk and surface exponents of the semi-infinite 
subsystems. 

These results have been tested through large scale 
Monte Carlo simulations, in which the critical behavior 
at the interface between 2D Ising, Potts and BW models 
was studied and satisfactory agreement has been found. 
However, it would be interesting to confirm the analyt- 
ical expressions for the interface exponents through a 
field-theoretical renormalization group study, using the 
methods of Rcf. 4. 

The results obtained in this paper can be general- 
ized into different directions. First we mention the case 
when the critical temperatures of the subsystems are 
not exactly equal but differ by an amount, AT C . If 
the deviation in temperature from the average value, 
T c = (T- + T+)/2, is small but satisfies T c - T > AT C , 
then our results are still valid. Our second remark con- 
cerns 3D systems in which sufficiently enhanced inter- 
face couplings may lead to an independent ordering of 
the interface above the bulk critical temperatures. In 
semi-infinite systems this phenomena is called the sur- 
face transition. 3 ^ 5 At the bulk critical temperature the 
ordered interface then shows a singularity, which is anal- 
ogous to the extraordinary transition in semi-infinite sys- 
tems. The singularities at the interface and extraordi- 
nary interface transitions remain to be determined, even 
in the mean-field approach. Third, we can mention that 
non-trivial interface critical behavior could be observed 
when one of the subsystems displays a first-order transi- 
tion. It is known for semi-infinite systems that the sur- 
face may undergo a continuous transition, which, how- 
ever, has an anisotropic scaling character, even if the 
bulk transition is discontinuous. 22,23 Similar phenomena 
can happen at an interface, too. Our final remark con- 
cerns the localization-dclocalization transition of the in- 
terface provided an external ordering field is applied. For 
two subsystems having the same ip 4 mean-field theory 
and the same 24 or different 25 critical temperatures, this 
wetting problem has already been solved. This solution 



could be generalized for subsystems having different field- 
theoretical descriptions. 
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APPENDIX A: SURFACE CRITICAL BEHAVIOR 

The calculation of the surface behavior is straight- 
forward when £/£ — > at the critical point. Here we 
give some details about the two cases where the surface 
boundary condition leads to a constant value for Z/£. 

1. ip 3 model with A > 

Using the results of Eq. (17), the boundary condition 
in (22) is rewritten as 



2A 



Let 



=f s K4)» sh i4)' (A1) 



3tanh 2 (^-)-l=at" 2 , #(0) = ft" 2 ■ (A2) 



cosh 



I 



(A3) 



Then, to leading order 

tanh (i) = 7I' -- V2C 

The first relation gives I in Eq. (25), and (Al) leads to 

2CA 



(A4) 



Finally, combining Eqs. (A2) and (A4), one obtains the 
value of ip(0) given in Eq. (26). 



2. ip G model with A < 

Since the surface is ordered, the profile is given by 
Eq. (17). The boundary condition in (22) translates into 



V2C S 
|A| 



Staiilr | -4 

2£ 



1 



x cosh — 
2£ 



3 tanh 2 ( 



1/2 3V2C / Z 
= sinn — 

3/2 

I . (A5) 
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The boundary condition is satisfied when 

J 3 tanh 2 (J^j - 1 = at 1 / 4 , (A6) 



so that: 



#(0 ) = ^ t -V4. 



(A7) 



Eq. (A6) gives the value of I in Eq. (31). Using the val- 
ues given in (A3) which remain valid here together with 
Eq. (A6) in (A5), one obtains: 



a _( 2C|A| 



1/2 



(A8) 



Inserting this expression in (A7) leads to the surface or- 
der parameter given in Eq. (32). 



Introducing this result in the second equation, one ob- 
tains an equation of the second degree in x = \fQ: 



. x ^ +2 m,m x+ ^=o. ( B4) 



Thus we have 



+ = / 



c_|A| to + 



This last result together with Eqs. (B3) and (B2) leads 
to the expressions given in (39) and (41). 



2. -»■ const., <C 1 



APPENDIX B: INTERFACE CRITICAL 
BEHAVIOR 

In this Appendix we give some details about the cal- 
culations of l± and p(0) limiting ourselves to three rep- 
resentative cases. Other results are easily obtained using 
similar methods. 



1. l±/£± <£ 1 

This situation is encountered for the p 3 — p 4 interface 
with A < and A — > oo as well as for the ip 4 — p 6 and 
p 3 — p 6 interfaces with A > and A — > oo. Here we 
consider as an example the p 3 — p 4 interface with A < 0. 

The boundary conditions in Eq. (38) are satisfied with 
ip-(z) and ip + (z) given by (16) and (18) and reads 



<^(0)=^ 1 / 2 coth(A- 



2 



3 cotlr 



l_ 
2~ 



-1 



C+p^t 1 / 2 . ZC_p t , (l_ 



: H—) 



Qp + t^ 
IAI 



coth 



2£+ 



(Bl) 



With l±/£± <C 1, one may expand the hyperbolic func- 
tions in powers of i±/(2£ ± ) To leading order, the first 
equation in (Bl) gives 



p(0) = 



3^* C% 
L 



so that 



l_ 



' 3ip l + t 



(B2) 



(B3) 



This behavior is obtained only for the ip 3 — <p 4 interface 
with A > 0. When A is smaller than a critical value 
A c to be determined later, the profiles <p>-(z) and <p+(z) 
are given by (17) and (19). They lead to the following 
boundary conditions: 



^ (0 )=^V2 tanh ( M = ^ 



2& 



3 tanh^ 



2£- 



■1 



With 



^tanhfA 



l_ 

2"e 

(B6) 



the first equation in (B6) gives 

p(0) = p+t = ^(3 tanh 2 K< 1) . 

It follows that 

l + = f®- (3 tanh 2 K < - l)£t . 
The second equation in (B6) can be rewritten as 

£ cosh K < A 



(B7) 



(B8) 



(B9) 



(BIO) 



Close to the critical point, the second term can be ne- 
glected so that 



(BID 
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Combining this result with (B9), one obtains 
'1 + 2A/A C A _ n Qp 



tanh_ft' < = 



2^o + - (B12) 



Since tanh if,, < 1, this solution remains acceptable as 
long as A < A c . Eqs. (B7), (B8), (Bll) and (B12) im- 
mediately lead to the expressions given in (44) and (46). 

When A = A c , l_ diverges and the order parameter 
remains constant, keeping its bulk value on the p 3 side 
of the interface. 

When A > A c , the profile is always increasing. Then 
<P-(z) is given by Eq. (16) and the boundary conditions 
are changed into 





_lp t 


3 coth 2 






2 







^1 cosh - 2 (I± 

2C+ V2^ 



3C_ p 1 



cosh( — — IsinlT 



^!tanh(^ 
A \2t + 



Inserting 



2£- 



K > 



(B13) 



(B14) 



into the first equation of (B6) leads to 

V(0) = ¥>J* (J^j = ^(Scoth 2 ^ - 1) (B15) 



and 



1+ = (3 coth 2 K > - 1)£+ 
fa 

From the second equation in (B13) one deduces 



(B16) 



acosh^ 1/2= r ^ (B17) 



£ sinh A> 



where the second term can be neglected close to the crit- 
ical point. Thus l + is still given by 



Comparing with (B16), one obtains 
coth K > 



1 + 2A/A C 



(B18) 



(B19) 



with the value of A c given in Eq. (B12). Since coth if-, > 
1, this new solution replaces the preceding one when 
A > A c . The results given in (45) and (46) follow from 
Eqs. (B14), (B15), (B18) and (B19). 



3. 1-/S,- -C 1, -»■ const. 

This is the situation encountered for the p 4 — p 6 and 
the ip 3 — (p 6 interfaces with A < 0. The treatment is simi- 
lar in both cases but we give some details for the p 3 — p 6 
interface which is a little more complicated. The inter- 
face is more ordered than the bulk so that the profiles 
are given by (16) for z < and (20) for z > 0. They lead 
to the following boundary conditions: 



p(o) = Vzp+t 1 ' 4 

3 coth 
sinh 



2 

Zy^C+tp+t 1 / 4 



3 tanh 
l_ 



2£ + 



+ 



1±. 

2£ + 

2£ + 
cosh 



cosh 



3 tanh 
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|A| 



1±. 

2£ + 

- 1 



2£+ 

1 -3/2 



-1/2 



1 



l_ 

3 tanh 2 



sinh 



2£+ 



l_ 

w 

- 1 



-1/2 

(B20) 



As in Appendix A 2, the solution is obtained by assuming 
that close to the critical point: 



3 tanh 



2£ + 



l = at 1 ' 4 . (B21) 



From this expression one deduces the leading contribu- 
tion to l + given in (57). 

With <C 1 the first equation in (B20) can be 

rewritten as 

P(0) = ^=6p-J^Y . (B22) 



Thus we have 



I- = Co 



(B23) 



The second equation in (B20) allows us to determine the 
value of a. Actually, we obtain the following equation 
for x = 1/y/a: 



x 4 + ax — b = . 



= 2 3 / 4 




^, 6=V3 



C + |A|- 



(B24) 



It is easy to verify that this equation has a single real 
positive root, xq . Below, we evaluate xq in the two 
limiting cases, x -C 1 and x ^> 1. 

When xq <C 1, one can iterate the relation 



x 



- 1-^ 



(B25) 
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following from Eq. (B24). Wc obtain 



priate when xq 3> 1 is 



x = 



1 + 



b 3 



23/4a|A|V^o + 



-1/2 



(B26) 



This result is valid as long as b 3 <§C a 4 — i.e., when 

i A '» A * = c '(S 2/a (|)' /S (l) 4/a - (B27 » 

Combining (B22), (B23) and (B26) one easily obtains the 
results given in Eqs. (57) and (59) for |A| > A*. 

The relation following from Eq. (B24) which is appro- 



x 
so that 

Xq 



1/4 



~ b 1 / 4 (l 



& V4 



"46" 



/o C*£o 
V CV|A| 



1/4 



a 



(B28) 



-1/2 



(B29) 



The correction term can be neglected when a 4 <C b 3 - 
i.e., when |A| < A*. Then Eqs. (B22), (B23) together 
with (B29) lead to the results given in (57) and (59) for 
|A| « A*. 
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